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Abstract. Optical emission-line profiles were evaluated in order to explore the structure of galactic nuclei containing 
H2O megamaser sources. Long-slit spectra of IC 2560, NGC 1386, NGC 1052 and Mrk 1210 were obtained at 
~ lOOkm/s spectral and ~ 2" x 2" spatial resolution. The following individual properties of the objects were 
found: The active nucleus of IC2560 (innermost ±2") emits lines typical for a high-ionization Seyfert-2 spectrum 
albeit with comparatively narrow profiles (FWHM ~ 200km/s). Line wings are stronger on the blue side than 
on the red side suggesting outflow. The central velocity gradient fits into the general velocity curve of the galaxy. 
Attributing it to a rotating disk coplanar with the galaxy leads to a Keplerian mass of ~ 10 7 M© inside a radius of 
100 pc. — The central 6" sized structure seen on HST Ha and [Om] images of NGC 1386 appears to be the inner 
part of a near-edge-on warped rotating spiral disk that is traced in Ha within a diameter of 17" along p. a. 23°. This 
interpretation is based on observed kinematic continuity and a typical S-shaped dust lane crossing the kinematical 
center. The central velocity gradient yields a Keplerian mass estimate of ~10 8 Mq inside R — 100 pc and 5 10 9 Mq 
inside 0.8 kpc. The total mass of the ionized gas of 10 5-6 Mq is small compared to the dynamical mass of the spiral 
disk. — The kinematical gradient of the rotating gas disk in the center of the elliptical galaxy NGC 1052 yields a 
mass of ~6 10 s Mq inside 166 pc. Two components are found: component A arises in a rotating disk, component B 
is blueshifted by ~ 400km/s relative to the disk. B likely originates from outflowing gas distributed within a wide 
cone. There is no need for a broad Ha component in unpolarized flux. The peculiar polarization seen in [Ol]A6300 
by Barth et al. (1999) may be related to the outflow component. — In Mrk 1210, a redshifted component 1Z 
contributes to the exceptional width of [Ol]A6300. An additional blueshifted component B strong in [Om] allows 
to fit the Ha + [Nil] blend without a broad component of Ha. B and 1Z are likely to be outflow components. The 
location of the brightness maximum BM is slightly (~ 1") shifted relative to the kinematical center KG Locating 
the true obscured nucleus in KC makes it understandable that BM displays faint broad scattered Ha in polarized 
light. — Galactic rotation and outflow of narrow-line gas are common features of this sample of water-megamaser 
galaxies. All decomposed line-systems exhibit AGN typical line ratios. Recent detections of H2O megamasers in 
starburst galaxies and the apparent asssociation of one megamaser with a Seyfert 1 AGN suggest that megamasers 
can possibly be triggered by optically detectable outflows. The frequently encountered edge-on geometry favoring 
large molecular column densities appears to be verified for NGC 1386 and IC 2560. For NGC 1052 and Mrk 1210, 
maser emission triggered by the optically detected outflow components cannot be ruled out. 

Key words. Galaxies: active - Galaxies: general - Galaxies: nuclei - Galaxies: starburst - Galaxies: individual: 
NGC 1052, NGC 1386, IC 2560, Mrk 1210 



1. Introduction 

During the past two decades, powerful H2O masers 
have been detected in the 6i6~523 transition at 22 GHz 
(A~1.3cm) towards almost thirty galaxies known to con- 
Send offprint requests to: H. Schulz, hschulz@astro.ruhr-uni- 
bochum.de 

* Based on observations collected at the La Silla outstation 
of ESO. The data have been evaluated with the MIDAS data 
analysis system (version Nov. 96) provided by ESO 



tain an active galactic nucleus (AGN) of Seyfert 2 or 
LINER type (dos Santos & Lepine 1979; Gardner & 
Whiteoak 1982; Claussen et al. 1984; Henkel et al. 1984, 
2002; Haschick & Baan 1985; Koekemoer et al. 1995; 
Braatz et al. 1994, 1996, 1997; Greenhill et al. 1997a, 2002, 
2003; Hagiwara et al. 1997; Falcke et al. 2000). These H 2 
masers (hereafter 'megamasers') appear to be more lu- 
minous, by orders of magnitude, then the brightest H2O 
masers observed in galactic star forming regions. All in- 
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tcrfcrometric data of H2O megamasers point toward a nu- 
clear origin. The emission originates from the innermost 
parsec(s) of the parent galaxy (e.g. Claussen & Lo 1986; 
Greenhill et al. 1995b, 1996, 1997b; Miyoshi et al. 1995; 
Claussen et al. 1998; Trotter et al. 1998; Herrnstein et al. 
1999). Adopting the so-called unified model (Antonucci & 
Miller 1985; Antonucci 1993) with Seyfert 2 nuclear disks 
being viewed edge-on, megamaser activity must be related 
to large column densities along the line of sight, a picture 
in line with the presence of dense molecular gas. 

There exist at least two different classes of mega- 
masers, those forming a circumnuclear disk (like in 
NGC4258) and those associated with nuclear jets (e.g. 
NGC1052). The observed H2O transition traces a dense 
(J5l0 7 cm -3 ) warm (Tki n Jb400K) molecular medium. The 
additional presence of jets and a nuclear X-ray source re- 
quires, however, the coexistence of neutral and ionized gas 
in the nuclear regions of these galaxies. In an attempt to 
study connections between megamaser emission on mil- 
liarcsecond scales and optical phenomena on arcsecond 
scales, here we present and analyze optical emission line 
spectra from four megamaser galaxies. 

2. Data acquisition 

The spectra were obtained in January 1996 using the 
Boiler & Chivens spectrograph attached to the Cassegrain 
focus of the ESO 1.52 m spectrographic telescope. A log 
of the observations, in particular the observed wavelength 
ranges, are given in Table 1. The detector was La Silla 
CCD No. 24 (FA2048L CCD chip with 15 //m wide square 
pixels). The spatial resolution element is 0'.'68 px -1 . Seeing 
and telescope properties limited the spatial resolution to 
the range 1'.'2 - 2" 5 which was determined by the width 
of standard star spectra on the focal exposures. The slit 
was visually centered and guided on the nuclear brightness 
maximum or a nearby offset star. The 2" wide slit projects 
to a spectral resolution of ~ 2A (~ 100 km s _1 ) as is con- 
firmed by the full width at half-maximum (FWHM) of 
comparison lines and the [Oi]A6300 night sky line. 

Employing the ESO MIDAS standard software (ver- 
sion Nov. 96) the CCD frames were bias subtracted and 
divided by well exposed flat-field frames. Night-sky spec- 
tra 'below' and 'above' any notable galaxy emission were 
interpolated in the region of the galactic spectrum and 
subtracted in each case. From such cleaned frames either 
single rows (corresponding to 0'.'68 and 2" along and per- 
pendicular to the slit direction, respectively) or the aver- 
ages of two adjacent spectral rows (1'.'36 x 2" formal spatial 
resolution) were extracted. 

Tug's (1977) curve was used to correct for atmospheric 
extinction. The spectra were flux calibrated using the 
standard star 40 Eri B (Oke 1974), but we use the cal- 
ibration only for relative line ratios due to occasional thin 
cirrus. Forbidden- line wavelengths were taken from Bowen 
(1960). Spectra in the blue and red region of the bright el- 
liptical NGC 3115 were taken as well and we used them as 
templates to subtract the stellar contribution from those 



spectra in which we analyzed the line-profile structure and 
intensity ratios. 

3. General data and morphology 

Table 2 (column 2) exhibits the morphological diver- 
sity of the four galaxies although we only refer to the 
crude Hubble type. While NGC 1052 is an elliptical galaxy 
(Fosbury et al. 1978), NGC 1386 (Sa according to Weaver 
et al. 1991) and IC2560 (Sb in Fairall 1986; SBb in 
Tsvetanov & Petrosian 1995) are disk systems whereas 
Mrkl210 is hard to classify from ground-based data. It 
looks roundish with an amorphous body supplemented 
by a few shell-like features, possibly a late-stage merger 
(Heisler & Vader 1994). An HST-snapshot image unveils 
a kpc-sized central spiral apparently seen close to face-on 
that can be classified as Sa in a Hubble scheme (Malkan 
et al. 1998), although it should be cautioned that this spi- 
ral may be gaseous rather than a density-wave triggered 
stellar spiral. 

Column 3 gives the heliocentric value of the systemic 
velocity as derived from our data. All velocity values given 
below are heliocentric. Column 4 of Table || gives distances 
(see footnotes) in accordance with recent work to ease 
comparison. Column 5 gives the linear scale in the galaxy 
applicable along the major axis for the distance (col. 4) as- 
sumed. Cols. 6 and 7 give the position angle of the appar- 
ent major axis and the angle of inclination from face-on. 
For the present work, these numbers should preferably cor- 
respond to the disk in which emission lines are generated, 
which is, however, only poorly determined. Comments and 
footnotes explain the origin of the numbers given. 

4. Results on individual galaxies 

Details of the spectral data analysis to find clues for basic 
kincmatical bulk components in the shape of the emis- 
sion lines are outlined in Appendices A - B (treatment 
of individual spectra) and Appendix C (determination 
of rotation curves). Here we discuss our results, starting 
with the two strongly inclined spiral galaxies IC 2560 and 
NGC 1386. Then we turn to the apparently more complex 
spectra of NGC 1052 and Mrk 1210. 

4.1. The Seyfert I C 2560 
4.1.1. Kinematics 

Along p. a. (90°, 270°) we traced regions with emission 
lines 22" east until 50" west. The reference position 0" on 
the abscissa of Fig. 1 marks the spot where the maxi- 
mal intensity of nuclear line emission is reached. Here 
^hei = (2922 ± 20)km/s is measured, which we identify 
with the systemic velocity. Outside the innermost nuclear 
region (useful line measurements are possible from the nu- 
clear maximum within a diameter of 5") appears a gap 
with no traces of a line and emission lines further out 
originate in disk HII regions. The corresponding velocity 
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Table 1. Journal of the spectroscopic observations 



Object/helioc. corr 




Frame 


Date (beg.) 


p. a. (slit) 


Pos. of slit 


A range (A) 


Exp. time (s) 


NGC 1052 


F020 


ll-Jan-96 


90° 


through nucleus 


5400-7400 


3600 


-27 km/s 


F056 


12-Jan-96 


90° 


through nucleus 


5400-7400 


1800 




F057 


12-Jan-96 


90° 


2" south of nucleus 


5400-7400 


1800 




F058 


12-Jan-96 


90° 


2" north of nucleus 


5400-7400 


1800 




F078 


13-Jan-96 


45° 


through nucleus 


5000-6800 


3600 




F107 


14-Jan-96 


45° 


through nucleus 


4000-6000 


3600 




F108 


14-Jan-96 


45° 


2" northwest of nucleus 


4000-6000 


3600 


NGC 1368 


F023 


ll-Jan-96 


90° 


through nucleus 


5400-7400 


3600 


-17 km/s 


F026 


ll-Jan-96 


90° 


2" north of nucleus 


5400-7400 


3600 




F080 


13-Jan-96 


23° 


through nucleus 


5000-6800 


2700 




F112 


14-Jan-96 


23° 


through nucleus 


4000-6000 


3600 


Mrkl210 


F066 


12-Jan-96 


90° 


through nucleus 


5400-7400 


3600 


+6 km/s 


F084 


13-Jan-96 


90° 


through nucleus 


5000-6800 


1333 




F085 


13-Jan-96 


90° 


2" south of nucleus 


5000-6800 


3600 




F087 


13-Jan-96 


90° 


2" north of nucleus 


5000-6800 


3600 




F116 


14-Jan-96 


90° 


through nucleus 


4000-6000 


3600 


IC 2560 


F030 


ll-Jan-96 


90° 


through nucleus 


5400-7400 


3600 


+20 km/s 


F118 


14-Jan-96 


90° 


through nucleus 


4000-6000 


3600 



a This heliocentric correction was added to the measured radial velocities. 



Table 2. Adopted basic data of the four galaxies (whei from this work) 



Object 


Type 


Vhcl 


dist. 


lin. sc. 


maj. ax. 


inclin. 


Comment 


Coord. (1950) 


Hubble 


km s _1 


Mpc 


pc/" 


o 







NGC 1052 


E 


1462 ± 20 


17.10 


82.9 


4c|j 


~ 6C0 


p. a. (major axis) and inclination 


0238-08 














of nuclear gas disk 


NGC 1386 


Sa 


878 ± 20 


19.80 


96.0 


23 


74 


major axis and inclination of galaxy 


0334-36 














on large scale (Tully 1988) 


Mrk 1210 


amorph 


4022 ± 30 


53.60 


260.0 


^ HO0 


« 


maj. ax. from gas kinematics; inclin. from image 


0801+05 












of kiloparsec spiral (Malkan et al. 1998) 


IC 2560 


SBb 


2922 ± 200 


260 


126.1 


45 


630 


p. a. and inclination of galaxy; distance 


1014-33 














39 Mpc if peculiar motion is ignored 



a distance from tWcss = 1282 km/s relative to CMBR (de Vaucouleurs et al. 1991) as used in Claussen et al. (1998) 

b Davies & Illingworth (1986) 

c from this work 

d we assume the same distance as for NGC 1365 in Fornax cluster; cf. Schulz et al. (1999) 

e distance via direct application of Ho = 75 k m s" 1 Mpc -1 

' kinematical major axis, this work; see Sect. 4.4.2 

9 type of inner kpc-spiral; based on HST image (Malkan et al. 1998) 

h Vhd measured by us at the location of maximal brightness; the value agrees with that of HI (Mathewson et al. 1992) 

! member Antlia cluster; distance from Aaronson et al. (1989) as used in Ishihara et al. (2001) 



major axis and inclination of galaxy on large scale (Mathewson et al. 1992) 



curve is shown in Fig.|l|. Between 15" and 20" west an 
interesting local up- and downturn (presumably density- 
wave related) can be noted in the spiral arm. This and 
two other intersections of a spiral arm (~ 45" to 50" west 
and 17" to 22" east) exhibit contiguous Ha emission and 
mostly [Nil] as well, with each region a few hundred parsec 
in extent. 

To evaluate the central linear gradient a deprojection 
onto the plane of the galaxy is required. Here 9 = 45° is 
the angle in the plane of the sky between slit position and 



major axis, because the slit extends along p. a. 90° and 
the galaxy's major axis along 45°, while its inclination 
i = 63° (Table |). 

Altogether, a nuclear gradient in rotational velocity of 
20km/s/(100 pc) is obtained yielding a Keplerian mass 
estimate (M = Rv 2 /G) of M loa = 9.310 6 M Q within a 
nuclear radius R — 100 pc. This is only 2.6 times the 
value that Ishihara et al. (2001) obtain for the mass inside 
the central 0.07 pc (note that they assumed an edge-on 
disk). Extrapolating our gradient to this radius yields a 
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negligible velocity compared to the actual velocities of the 
maser sources (210-420 km/s) measured by Ishihara et al. 
(2001) between 0.07 and 0.26 pc. This might be explained 
by an inner Keplerian upturn of the rotation curve due to 
the central mass concentration. 

The line profiles from the nuclear region within ±2" of 
the center (defined by the maximum of continuum and 
line brightness) are relatively narrow for a Seyfert, with 
FWHMs of (200-250) km/s. Fig. | shows that the blue side 
of the [Oiii] profiles can be fit by a Lorentzian while Fig. [| 
shows that the red side of the [Om] profiles can be fit by a 
Gaussian (this also holds for H/3, Ha and [Nil]). This indi- 
cates that the observed [Oiii] profiles have a stronger outer 
wing on the blue side, which is most likely a signature for 
outflow. However, the amount and velocity of the outflow 
cannot be modelled precisely because the line width is only 
twice our spectral resolution. Taking the Gaussian as rep- 
resentative for the rotational component convolved with 
the instrumental profile, the deconvolved l.o.s. outflow- 
velocity component is estimated to be ~ 100 km/s. 

It should be noted that our near-edge-on view of the 
central kpc-disk may be unfavorable to measure outflow- 
ing gas, because outflow is frequently directed perpendic- 
ular to a galactic disk. 



Velocity curve of IC 2560 




5020 5040 
X/ANGSTR. 



Fig. 2. [Om]AA4959, 5007 from the central 1"36 (E-W) x 
2"(N-S) of IC 2560 (F118) (plotted in the observed frame 
as all subsequent spectra). The Lorentzians used for the 
fits are too wide within the red wings. 




5020 5040 
VAngstr. 



Fig. 3. The same [Oiii]AA4959, 5007 observations as in 
Fig. 2, but now fit by Gaussians, which turn out to be 
too narrow to fit the blue wings. 



Fig. 1. Heliocentric radial velocities (averages of the peaks 
of Ha and [Nn]A6583) measured in spectrum F030 of 
IC2560 (along p. a. 90°; abscissa in arcsec). Statistical un- 
certainty ±11 km/s. 



4.1.2. Trends in line-ratios 

For readers interested in line-excitation variations a few 
trends of line ratios are given in the present, albeit primar- 
ily kinematic, study. In the innermost ±2" along p. a. (90°, 
270°) the intensity ratios of [Nil] /Ha and [Sn]/Ha are in- 
creasing to the west. At the three positions, central maxi- 
mum, 1" west, and 2" west, we measure: [Oi]A6300/Ha = 
0.15, 0.2 and 'not measurable'; [Nil] A6583/Ha = 1.05, 
1.2 and 1.5; [Sn](A6716 -I- A6731)/Ha = 0.6, 0.65 and 
0.9, respectively. [Om] was extracted over 2" bins yield- 



ing [Oin]A5007/H/3 = 13 in the central bin, which de- 
creases to 10.5 at 2" west, while it may slightly increase 
to 14 at 2" east. The other line ratios at 1" east and 
2" east do not significantly differ from the central val- 
ues. The density-indicator [Sn]A6716/[Sll]A6731 is found 
in the range 0.90 ± 0.05 in this inner region, correspond- 
ing to single-component densities in the range (7— 10) 10 2 
cm~ 3 . 

The relative strengthening of the low-ionization lines 
to the west (which is also the far side if the spiral is trail- 
ing) may be due to a decrease of the ionization parameter 
by increased filtering of the ionizing radiation with nu- 
clear distance. Such an effect is demonstrated in Schulz 
& Fritsch (1994). The directly measured line intensities of 
Ha have a steeper fall-off to the west than to the east. 

Exploratory spectra with an east-west oriented slit po- 
sitioned 2", 4" and 6" south of the nucleus (taken to look 
for some circumnuclear influence of the activity) are too 
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weak for detailed measurements. At 4" south the spec- 
trum appears to be of lower excitation than known from 
the standard LINER sample shown in the diagnostic dia- 
grams of Schulz & Fritsch (1994). 

4.2. The Seyfert NGC1386 

4.2.1. Global kinematics of the inner emission-line 
region 

We here first describe a velocity curve derived from 
the line peaks of Ha and [Nn]A6583 measured along 
the large-scale major axis at p. a. 23° (Fig. Q). The 
[Oiii]AA4959, 5007 (spectrum F112) peak positions agree 
in the innermost ±3" where they can be measured. For a 
comparison o f the complicated line profiles of [Oiii] and 
Ha see Sect. 1.2. 2|. H /3 is faint, partially due to strong 
extinction (Sect. 4.2.2 ), and cannot be measured beyond 
~ 2'.'5 of the brightness maximum. 



n i i i i r 



NGC 1386 p. a. 23 



600 1 — L 
-15 



-5 5 

ARCSEC 



Fig. 4. Velocity curve of NGC 1386 along p.a. 23°(from 
F080). The data correspond to the mean position of the 
line peaks of Ha and [Nn]A6583 in the raw spectrum. 

The zero point of the abscissa of Fig.[| corresponds 
to the spectral row where continuum and intensities of 
Ha and [Om] are at maximum. We call this brightness 
maximum 'region £>' because of its blueshiftQ relative to 
the kinematical center KC (see Fig. KC is defined by 
taking the symmetry center of the outer portions of the 
velocity curve. For its determination we took means of 
the flat portions of the velocity curve beyond ±8" yield- 
ing (after symmetrizing) v — 168km/s (northeast) and 
v = — 168km /s (so uthwest) around vkc = 877km/s. 
Below in Sect. 4.2.2 we will find that the mean of the 



line cores of six different lines arising within the inner- 
most 2" of KC yields 878km/s, which is in (fortuitously) 
perfect agreement with ukc- 

Linearly interpolating in the spatial direction, KC is 
found to lie — 0'.'6 (to the northeast) offset from B. 

1 Blueshifted regions in NGC 1052 and Mrk 1210 will be 
called B as well 



The position of B on the sky can be inferred from 
an HST Ha image shown in Ferruit et al. (2000) (their 
Fig. 4, left panel, middle row; hereafter FWM2000). The 
following description refers to their Ha image. Our spec- 
trograph's slit encompasses the bright parts of the image 
by a 2" wide band along p.a. 23°. The major axis of the 
Ha emission string appears to extend along p.a. 6°, or 
A9 = 17° offset from the slit position angle. 

From the rotation pattern and the appearance of 
NGC 1386 as a highly inclined galaxy one expects to look 
upon a similarly inclined emission-line disk in the cen- 
ter. Images of galaxies at high inclination commonly show 
a midplane-dust lane. In FWM2000's Ha image appears 
indeed a dark lane crossing the center between the two 
brightest spots that seems to wind from the southwest to 
the northeast in a extended S-shaped pattern. 

The dominating spot in FWM2000's image, just below 
the dust lane, has to coincide with B. B is located at the 
zero position in our Fig.^J. The kinematical center KC, 
offset by 0'.'6 in our spectrum, has to be located within 
the dust lane. 

The Ha image structure is characteristic for a near- 
to-edge-on and warped two-armed spiral disk. The spiral 
bending is in the same sense as that of the large galaxy, but 
the inclination appears to be larger than i — 74°, the in- 
clination of the large scale disk (Table 2). Since the north- 
eastern part is receding, the northwestern arm should be 
on the near side if we make the usual assumption that the 
arms are trailing. 

The above interpretation of the main extended struc- 
ture in FWM2000's Ha image may be considered to be 
at variance with the usual notion of a 'radiation cone 
along the minor axis'. When interpreting IR spectroscopy, 
Rcunanen et al. (2002) use the term 'radiation cone' to 
describe the morphology of the emission-line gas, but, on 
the other hand, they note that the H2 "rotation curve is 
relatively ordered parallel to the cone, suggesting that cir- 
cular motions dominate the kinematics" . Our data agree 
with this statement on the kinematics, but it is not clear 
whether the edge-on disk is ionized by a UV 'radiation 
cone' (small plumes along p.a. 90° could be due to out- 
flowing gas illuminated by a radiation cone; see end of 
Sect. 4.2.2). 

The interpretation in terms of dominating circular mo- 
tions of the inner ionized gas is further supported by kine- 
matical continuity with gas at larger galactocentric dis- 
tances. The bright Ha features seen between 4" northeast 
and 2" southwest in the HST image are kinematically 
smoothly joined with gas further out (that is traced in 
our spectra but not in the HST image) to 14" northeast 
and to 14" southwest (Fig. [!]). Furthermore, the map of 
peak velocities in Weaver et al. (1991; their Fig. 5) shows 
a slightly distorted, but nonetheless typical pattern of a 
rotating disk. An interpretation by a bipolar flow acci- 
dentally aligned with the kinematical major axis of the 
rotating disk would appear contrived, because it had to 
be essentially decoupled from the gravitational field and 
its accordingly stratified interstellar medium. 
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Altogether, the natural explanation of our Fig.[| is 
that of a slightly projected rotation curve, although the 
rotating rings may not have a common plane, because 
FWM2000 noted changes of cllipticity and p. a. (major 
axis) with galactocentric radius. Owing to its AGN-typical 
line ratios the Ha gas is ionized by nuclear AGN radia- 
tion, but since the disk is likely to be matter-bounded 
the UV radiation is not necessarily confined to a narrow 
cone. Therefore we dismiss the term 'radiation cone' for 
NGC 1386. 



For estimates of the rotation velocity utilizing Eq. C.l 



we have to know the orientation of the supposed single 
disk. Regarding apparent distances a projection correc- 
tion would be small (because cos 17° = 0.96). A velocity 
correction with i = 74° and 6> = 17° would formally yield 
a factor 1.5 to obtain a rotational velocity. However, this 
procedure is questionable because i = 74° corresponds to 
the outer galactic disk and the inner disk addressed here 
appears to be more inclined. 

For Keplerian mass estimates (that are good to a factor 
of two even if the potential is flattened; cf. Lequeux 1983) 
it is therefore sufficient to assume that the slit extends 
along the major axis and that the disk is seen edge-on, i.e. 
measured velocities are rotational velocities and angular 
distances along the slit correspond directly to linear dis- 
tances (an inclination correction solely with sin 74° would 
lead to 8% larger masses). 

The HST image (FWM2000's Fig. 4) only shows the in- 
ner 6" of the disk where Ha is brightest. As noted above, 
the velocity curve in our Fig.|] extends into the galactic 
disk beyond the bright Ha clouds seen in the HST image. 
The curve becomes flat at ±8'.' 5 corresponding to a galac- 
tocentric radius R = 816 pc at which v = 168km/s. These 
numbers yield a mass of M — Rv 2 /G — 5.4 10 9 M© in the 
central region of 1.6 kpc diameter of NGC 1386. In the 
center, the steep rise of the rotation curve to ~ 90km/s 
at 1" yields a mass of 1.8 10 s M© within a diameter of 200 
pc. 

From the observed Ha luminosity Ln a « 10 erg/s 
(Colbert et al. 1997) and adopting Lh q = 2.9 Lh/3 one 
can estimate the volume emission measure Vn 2 e for fully 
ionized (n e = n p , the proton density) hydrogen gas via 
Lh/3 = (Vrig) (471711,8 /«e) (jh/3 is the emission coefficient 
of H/3 per steradian). Taking from Osterbrock (1989; 
Table4.2) 4irj Hp /(n e n p ) = 1.24 10~ 25 erg cm 3 s" 1 (for 
T = 10 4 K) and adopting a typical NLR electron-density 
range n e — 10 2-3 cm -3 , applying a conservative extinc- 
tion correction of a factor of 3 and a factor of 1.45 for the 
presence of helium, one obtains 10 5-6 M© for the mass of 
line-emitting gas, which is small compared to the dynam- 
ical masses given above. 



4.2.2. Line profiles from the inner disk 

For sampling in closer accordance with the seeing we 
binned pairs of the 0'.'68 wide CCD rows to spectral 
scans, which each encompass light from 1'.'36 along the 



entrance slit, and subtracted scaled spectra of our tem- 
plate NGC 3115 to remove the contribution of stellar light. 
H/3 turned out to be faint and too noisy for an accurate 
estimate of the extinction from the Balmer decrement. 
Roughly, at the position of the line maximum we have 
I (Ha) / 1 (H(3) » 4, while this value approaches 6 at l'/36 
southwest, and 7 at l'/36 northeast. In a simple screen 
model these values translate into Ay = 0'. n 9, 2 rn and 2™5, 
respectively. 

Figs.|| and || show the 'blue' line profiles of [Oiii], and 
Figs, [f] and || display the 'red' line profiles of Ha and [Nil] 
from the 'central' (max. line emission) 1'.'36 along the slit 
and the next two l'/36 wide bins towards the southwest 
(p.a. (23°+ 180°)) and towards the northeast (p.a. 23°), 
respectively. Although the blue and red frames were ex- 
posed at l'/4 and 2'.'5 seeing, respectively, the profiles of 
[Oiii] on one side and Ha or [Nil] on the other side reveal 
the same basic shape. The most pronounced difference ap- 
pears at l'/36 southwest, the sharp shoulders on the red 
side of [Oiii] are smoothed out in Ha and [Nil] , which we 
trace back to the seeing effect. 




Fig. 5. [Oiii]A4959 and A5007 line profiles extracted from 
l'/36 wide spectral scans (of spectrum F112). Continuous 
line: from region B (line maximum); dashed: 1'.'36 south- 
west; dotted: 2'/72 southwest. 

In the last subsection we concluded that the spectrum 
along p.a. 23° arises in an inner gaseous disk that is seen 
near to edge-on and whose line emitting regions are en- 
compassed by the 2" wide slit. Hence, we expect line pro- 
files resulting from line-of-sight integrations through the 
disk. 

All profiles have pronounced narrow peaks (called line 
cores here) with steep flanks. The velocities of the narrow- 
line cores are given in Table||. Our spectral resolution 
is about half of the apparent width of the line core of 
disk-component B. The first numerical column in Table|| 
gives the velocities of region B, as measured with differ- 
ent lines. As Figs.|| to g or Tablc|| show, the radial ve- 
locity of the core of B is close to those at 1'.'36 south- 
west and 2'.' 72 southwest, while the regions to the north- 
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Fig. 6. As in Fig. || but now towards northeast. 
Continuous line: from region B (line maximum); dashed: 
1'.'36 northeast; dotted: 2'.'72 northeast 



NGC 1386, Ha, [Nil] to NE 




6540 6560 6580 6600 6620 

X/Angstr. 

Fig. 8. As in Fig. [7], but now towards northeast. 
Continuous line: from region B (line maximum); dashed: 
1'.'36 northeast; dotted: 2'.'72 northeast. 



NGC 1386, Ha, [Nil] to SW 1 




Fig. 7. Ha+[Nn] AA6548, 6583 line profiles extracted from 
l'/36 wide spectral scans (of spectrum F080). Continuous 
line: from region B (line maximum); dashed: l'/36 south- 
west; dotted: 2'.' 72 southwest. 



east are redshifted. Hence, region B is on the 'blue' - 
approaching - side of the rotating disk. In Figs.|| and 
H we see the core from B and the cores from the northeast 
line profiles together. The radial velocity of the center of 
the disk (the disk-center velocity) should then be amidst 
the southwest and northeast cores, at u£ el = 878km/s 
(mean of six spectral lines). This is in agreement with 
vkc = 877km/s determined from folding the outer rota- 
tion curve (Sect. 4.2.1). 

Schulz et al. (1995) (their Fig. 4) showed that a ro- 
tating edge-on disk with an emissivity law falling with 
radius will produce asymmetric profiles, with the more ex- 
tended wing pointing towards smaller absolute values of 
rotation velocity. However, the wings should not pass the 
disk-center velocity unless there are strong random mo- 
tions and/or radiation 'from the other side' is scattered in 
via observational seeing (and guiding inaccuracies), as is 
shown by an example in Sect. 7.2 of Schulz et al. (1995). 



The asymmetric line profiles of component B have the 
extended wing on the side expected for a rotating disk, 
but their shoulder extends to longer wavelengths than the- 
oretically predicted unless there is strong spatial smear- 
ing or an inner high-velocity branch along the same l.o.s. 
(e.g. central Keplerian rotation). The shoulders even ex- 
tend beyond the northeast-cores (Figs. |]and ||). Note that 
in velocity units the red flank of the core l'/36 northeast 
lies at about the upper limit of the rotation curve in Fig.|| 
(~ 1050 km/s) and the blue flank of the core of compo- 
nent B (or the core of the component l'/36 southwest) 
lies at the lower limit of the rotation curve (~ 700 km/s). 
Consequently, the outer red shoulder of profile component 
B cannot be explained by light from regions of the outer 
rotation curve, while the sharp transitions to the inner 
part of the red shoulder (particularly in Ha and [Nil] of 
B) and in [Om] 1'.'36 southwest suggest that this is emis- 
sion from the northeast side of the disk. This picture im- 
plies that there is nearly no emission seen from the nucleus 
of the disk at the disk-center velocity (in agreement with 
the location of KC within the dust lane; see Sect. 4.2.1 ) 
as such emission would have smoothed out the transition 
from the southwest to the northeast component. 

Outer wings are found on the blue and red side of com- 
ponent B and extend up to 900 km/s symmetrically to the 
disk-center velocity in [Om] , but not symmetrically to the 
core velocity of B. Wings can also be traced to a slightly 
lesser extent from the components at 1'.'36 and 2'.'72 dis- 
tance from B. In [Nil] outer wings appear to extend to 
- 700 km/s. 

A possible scattering origin of such wings will be dis- 
cussed in Sect. 5.1.2. Since the kincmatical center is ob- 



scured by the dust lane, an explanation by fast Keplerian 
motion close to the disk center is unlikely. 

Spectrum F023 taken along p. a. 90° (along which 
FWM2000 detected two emission-line plumes) shows emis- 
sion over ±2" with a slightly declining velocity gradi- 
ent from east to west (~ 5 km/s/"). The near-north-south 
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Table 3. Measured velocities (km/s) of line cores along p. a. 23° in the center of NGC 1386. Uncertainty 15km/s 



Line At max. (B) l'.'36 southwest 2'.'72 southwest l"36 northeast 2'.'72 northeast 



[Om]A4959 776 757 747 964 983 

[Om]A5007 793 762 746 962 984 

H/3 777 767 751: 965 965: 

[Oi]A6300 849: 750 741 947 996 

Ha 790 777 770 967 993 

[Nn]A6583 789 775 768 980 1003 

[Sn]AA6716,6731 798 776 771 979 1016 



orientation of the disk explains why a notable east-west 
gradient in the l.o.s. projections of rotational motions is 
unexpected. Its near-edge-on orientation makes the detec- 
tion of minor-axis outflow difficult. The gradient could be 
due to outflow along the plumes if the disk is slightly tilted 
from edge-on. The extended red wing in Ha, [Nil] and [Sn] 
(seen in the F080 spectrum of component B (continuous 
line) in Figs.^ and ||) is confirmed by F023. 

4.2.3. Line ratios 

Albeit of the primarily kinematic nature of the present 
study, a few relative intensities of the major diagnostic 
lines arc to be noticed in order to identify the nature of 
the emission-line region in question. Weaver et al. (1991) 
noticed a wealth of varying emission-line ratios from the 
typical Sy-2 high-excitation ratios to extranuclear LINER- 
like ratios. Even in the few strips that our spectra cover, 
drastic line-ratio changes are seen. Along the p. a. (23°, 
23°+180°) slit, [Nn]A6583/Ha, [Sn](A6716+A6731)/Ha 
and [Oi] A6300/Ha increase from 1.5, 0.68 and 0.13 at the 
brightness maximum (B) to 2.5, 1.10, 0.23, respectively, 
at 4" southwest. To the northeast there is only a moder- 
ate (10%) increase towards 3" northeast in [Nil] and [Sn], 
only [Oi] rises to 0.2. 

[Olll]/H/3 is 13 at B, but rises to 21 at l'/36 southwest, 
and 17 at l'/36 northeast. [Sn] A6716/[Sn] A6731 is 1.0 at 
B (indicating n e — 7 10 2 cm~ 3 ) and 1'.'36 southwest, but 
rises to 1.2 for the next 3" on either side (indicating an 
average n e = 2.5 10 2 cm -3 ). 

A faint spectrum along p. a. 90° (F026), with the 
slit set 2" north of B, shows essentially constant values 
- 1 and - 0.8 (n e w 10 3 cm -3 ) for [Nn]A6583/Ha and 
[Sn]A6716/[Sn]A6731 over about 10", respectively. 



4.3. The LINER NGC 1052 

4.3.1. Fitting line profiles with single components 

Davies & Illingworth (1986) gave p. a. 49° as the approx- 
imate position angle of the major axis of the nuclear 
gaseous disk in NGC 1052. For deprojecting the observed 
velocities we assume that spectra F078 and F107 (slit 
along p. a. 45°) were taken along the major axis. 




Fig. 9. Heliocentric velocities/ (km/s) along p. a. 45°in the 
center of NGC 1052 (northeast to the left). The filled oc- 
tagons show the mean positions of the line peaks of Ha and 
[Nn]A6583 as measured in the raw spectrum. The open 
octagons display the mean values of single-component 
Lorentzian fits given in Table |^. 



Fig. [9] shows two major-axis velocity curves; one shows 
the positions of the peaks of Ha and [Nn]A6583 in the 
raw spectra from 0'.'68 row extractions. The other curve 
corresponds to template-corrected spectra binned to 1'.'36 
wide spectral extractions and displays the mean values of 
single- component Lorentzian fits (see Table^J) of the lines 
Ha, [Nll]AA6548, 6583, and [Sll]AA6716, 6731 from the disk 
(see Appendix B regarding the fitting procedures with 
Gaussians and Lorentzians as basic functions). Beyond the 
steep central gradient the curves flatten out. They are not 
symmetric relative to the zero point, which is located at 
the brightness maximum. 

Davies & Illingworth (1986) identified a nuclear ionized 
disk that emits most of the emission-line spectrum. For 
subsequent mass estimates an inclination of the disk i = 
60° is assumed as derived in the following subsection (note 
that Bertola et al. (1993) obtained i — 53° from the data 
of Davies & Illingworth; this value would increase the mass 
estimates given below by 18%). 

The central velocity gradient (between maximum and 
two adjacent pixels) of 85 km/s/" translates formally into 
59.2 km/s/(50pc), yielding a Keplerian mass of 4.1 10 7 M Q 
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Fig. 10. Red spectrum (including [Oi]AA6300, 6363, 
Ha+[Nn]AA6548,6583, and [Sn]AA6716, 6731) from the 
brightness maximum (1"36 x 2") of NGC 1052 with the 
slit aligned along p. a. 45°. Fits with single Lorentzians 
(dashed) or Gaussians (dotted) are insufficient. Positions 
and line widths from the fits with Lorentzians are given 
in Table ^ (column headed by 0) . 




6550 6600 
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Fig. 11. Two narrow-line systems A and B, each char- 
acterized by three Lorentzian functions with fixed line 
width and position in velocity space (the intensity ratio 
of A6583 to A6548 is fixed to 3 to 1), are sufficient to fit 
the observed Ha+[Nn]AA6548, 6583 blend from the region 
with maximum emission (continuous line shown above) 
and till 2" northeast. Ha+[Nn] from system B (dotted) 
and the 'fitting' sum A + B (dashed) are shown. While 
B is blueshifted and indicates outflow, A is consistent to 
arise in a rotating disk. 

inside a radius of 50 pc, which corresponds to the inner, 
just resolved core with 1'.'2 diameter. 

The processed data (Table |j) have 1'.'36 resolution and 
lead to a gradient of llOkm/s per 2" radius (measured at 
a velocity of 1478 km/s), which translates into 6.3 10 s M Q 
inside a radius of 166 pc. 

Albeit crude, these mass estimates imply that 
NGC 1052 has a core of relatively small mass. In order 
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Fig. 12. As in Fig.|Tl], but now for [Sn]AA6716, 6731. 
Component B is additionally constrained to be in the high- 
density limit (n e 3> 10 3 cm -3 ). 
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Fig. 13. As in Fig. |TT|, but now for [Oi]AA6300, 6363. 

to judge the reliability of these numbers a closer look on 
the fits is needed. 

The data given in Table | were obtained by the 
Lorentzian fits of [Oi]AA6300, 6363, Ha, [Nn]AA6548, 6583, 
and [Sll]AA6716,6731 shown in Fig.[lO| Although fitting 
the major part of the profiles these fits are unsatisfactory 
in several respects. 

Firstly, measured outer line wings appear too weak or 
too strong as compared to Lorentzian or Gaussian fits, 
respectively. In itself, this might not be taken as serious 
because it could simply reflect the deficiency of these func- 
tions in comparison to the typical core-wing structure of 
narrow-line profiles. 

Secondly, the final summed fit lies systematically below 
the bridges connecting [Nn]A6548 and Ha or [Nn]A6583 
and Ha, and connecting the [Sn] lines. It also lies below 
the blue wing of [Nn]A6548. This holds irrespective of the 
type of profile function. 

Thirdly, although no conspicuous asymmetry is seen in 
[Ol] , in this case the fit function turns out to be blueshifted 
by ~ 100 km/s and appears 80% wider with respect to 
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Fig. 14. As in Fig. [fl], but now for H/3 and 
[Oiii]AA4959,5007. 



other forbidden lines. Also, sometimes Ha requires con- 
siderably larger widths than [Nil] although its position is 
not significantly shifted (Table |]) . 

These findings indicate that the line profiles are com- 
plex and that a one component model is not sufficient 
for an appropriate simulation. In the next subsection we 
therefore decompose the line profiles in a more detailed 
manner. 

4.3.2. Decomposition of the line profiles 

In view of the problems encountered with single- 
component fits, a significant improvement may be ob- 
tained by a mixture of two line systems, conveniently de- 
noted as A and B, with system B blueshifted with respect 
to A. 

In Fig.|ll| a fit to Ha + [Nil] from the 'central' (center 
defined by maximum emission) 2" x 0'.'68 (slit width x 
width of CCD extraction row) is shown along p. a. 45°. 
The line blend can indeed be successfully decomposed into 
two narrow-line Ha + [Nil] systems. However, such a fit 
of a smeared out blend of lines can only be found by the 
software if guided by strong constraints, which we are now 
going to justify. 

In Table|B| final results for the FWHMs, velocities and 
line ratios of systems A and B are given. The data are 
restricted to the spatial regime between 2" northeast and 
0'.'68 southwest for which component B can be unambigu- 
ously recognized. 

The Ha+[Nn] blend has been fit by assuming that all 
three lines in each system are constrained to have the same 
radial velocity and line width w (=FWHM; see Appendix 
B) while the intensities of Ha and [Nil] are left free (the 
3:l-[Nn] ratio is fixed; see Appendix B). Each of the two 
systems is fixed by four parameters (same position and 
same width in velocity space, height of Ha and height of 
[Nn]A6548). Thus, only eight parameters are fit for six line 
profiles which would be determined by eighteen parame- 
ters in total if each single profile were free of constraints. 



A slight shoulder on the blue wing of [Sn]A6716 re- 
quires the presence of the blueshifted component in [Sn] 
as well. An example for a fit of the [Sn] system employing 
analogous precepts as for Ha + [Nil] (adopting w(A) and 
w(B) from the Ha + [Nil] fit of each region) is shown in 
Fig. [l2]. Five out of twelve parameters were fit (two ve- 
locities, for systems A and B, two amplitudes for system 
A and one amplitude for system B, where the [Sn] inten- 
sity ratio was fixed to 7(6731)//(6716) = 2.3, which is the 
saturation limit for electron densities n e ^> 10 3 cm -3 ). 

A large electron density of B has been chosen to in- 
crease the contribution of [Oi] AA6300, 6363 relative to 
[Sn]AA6716,6731. The critical density of the upper [Oi] 
levels, 1.7 10 6 cm -3 , is three orders of magnitude above 
that for [Sn] . A high electron density is indicated by the 
single-component fits which yield large relative strengths 
of [Oi] in the diagnostic diagrams of Veilleux & Osterbrock 
(1987; hereafter VO) and larger line widths than those 
of [Sn] and [Nil] (Table |). The trend that [Oi] can be 
boosted mainly by density is shown in Komossa & Schulz 
(1997; their Fig. 4, upper panel). Even though a distinct 
blue peak or a clear shoulder is missing in [Oi]A6300 
(Fig. 13), the two-component fit with Lorentzian functions 
is satisfactory and explains the greater width of the [Oi] 
profiles by a relatively strong B contribution. 

Due to the likewise high critical density of 7 10 5 cm -3 
of the 1 D 2 upper term of [Oin] AA4959, 5007, these 
lines may be density-boosted in high-ionization zones. 
Consequently component B might be conspicuous in these 
lines as well which is indeed verified by a significant blue 
shoulder of the profiles (Fig.|l4|). The [Oi] and [Oill] spec- 
tra were consequently fit in the same way as the [Sn] lines; 
making use of 3: 1-line ratios in each of the oxygen dou- 
blets, however, four free parameters out of a total of twelve 
were sufficient, one amplitude and one velocity for each of 
the two systems, A and £>, respectively. 

The row 'mean' in Table|| gives the average velocities 
over the four line groups (H/3+[0iii], [Oi], [Sn], Ha+[Nll], 
the latter is weighted twice) for velocity systems A and B. 
The line velocities for system A given in Table I are con- 
sistent with arising in a rotating circumnuclear disk (as 
suggested by Davies & Illingworth 1986). The differences 
relative to 1502 km/s (next row) exhibit a nearly linear 
gradient to the northeast of 57 km/s per 0V68 step (=56.37 
pc) or 51 km/s at a formal radius of 50 pc from the kine- 
matical center (which appears to be slightly shifted from 
the location of maximum emission). This would lead to a 
Keplerian mass estimate of 3 10 7 M inside 50 pc if the 
disk were edge-on, or 4 10 7 M for i = 60° as given below. 
The latter value is in agreement with the single compo- 
nent fits given in Sect. 4.3.1 because the blue component 
B hardly effects the overall fits (except for [Oi] which was 
not considered in the mean values of Table |). 

Two slightly poorer spectra along p. a. 90°(F020 and 
F056) were decomposed in a similar way (with outflow 
fixed by 400 km/s). From these a component-,4 gradient 
of 15 km/s per 0'.'68 was found. With p. a. (major axis) 
= 45° we can adjust the velocity deprojection factor of 



Hartmut Schulz & Christian Henkel: Rotation and outflow in water-megamaser galaxies 



11 



Table 4. This table gives heliocentric velocities and (after the |) values of w=FWHM (in km/s) derived by single- 
component Lorentzian fits in template-corrected 1'.'36 wide spectral extractions along p. a. 45° of NGC 1052. For [Oi] 
and [Nil] the doublet ratios 3:1 and the requirement of same widths (as in [Sn] ) additionally constrain the fits. The 
column labeled '0' refers to the 1'.'36 wide spectral scan at maximum emission. Last row: mean values of the lines Ha, 
[Nil] AA6548, 6583, and [Sn]AA6716, 6731 (open octagons in Fig.|). 



Line 


5'.'44NE 


4'.' 08 NE 


2'.'72NE) 


l'.'36NE 





l'/36 SW 


27 72 SW 


4'.'08 SW 


5'.'44 SW 


Hq 

[Nil] 
[Sii] 

[0 i]A6300 


1330|161 
1363 168 
1342|178 


1369| 137 
1381 128 
1374 156 
1361 360 


1369|226 
1366|205 
1374 229 
1325|401 


1376|492 
1380 316 
1391 365 
1330 647 


1430 [706 
1446 435 
1471|454 
1373|711 


1531|545 
1552 383 
1550 393 
1469 [709 


1594|270 
1607 299 
1605 264 
1622|331 


1615|196 
1631 229 
1645 208 


1614|172 
1629 201 
1621 200 


mean 


1345| 169 


1375| 140 


1370|220 


1382|391 


1449 [532 


1544|440 


1602|278 


1630|211 


1621|191 



Table 5. Line widths w=FWHM, heliocentric velocities and line ratios measured for the rotational component A and 
the outflow component B (separated by a |) along p. a. 45° from 0'.'68 southwest to 2'.'04 northeast in NGC 1052. 





2" 04 NE 


1736 NE 


0'.'68 NE 





0768 SW 


w{A)\w{B) 


268|208 


296|366 


378|435 


420 1 498 


340 1 565 


Line 


H/3, [Om] 


1354 931 


1407| 1059 


1468|1064 


1496|1081 


1629|1342 


Ha, [Nil] 


1371 998 


1414 1043 


1449 1066 


1505 1096 


1567)1249 


[Sii] 


1369 - 


1414 1109 


1448 1066 


1498|1154 


1559 1353 


[Oi] 


1355 992 


1397 1016 


1484 1141 


1506|1165 


1575|1256 


mean 


1364|980 


1409|1054 


1460| 1081 


1502|1116 


1579| 1290 


mean— 1502 


-138| - 522 


-93| - 448 


-42| - 421 


0| - 386 


77| - 212 


I(A5007)//(H/3) 


1.97|1.94 


1.74|2.69 


2.03|2.43 


1.71|2.82 


2.03|2.15 


7(A6583)//(Ha) 


1.16 1.37 


1.12|1.32 


1.05|1.09 


1.090.86 


1.23 0.89 


7([Sll] to t//(Ha) 


1.29 - 


1.130.90 


1.180.84 


1.36 0.62 


1.58 0.91 


I(A6300)//(Ha) 


0.48 1.79 


0.58 0.57 


0.480.87 


0.45 0.66 


0.39 0.59 


7(A6716)//(A6731)(.4) 


1.40 


1.15 


1.04 


1.17 





Eq. C.l so that along both slit directions the rotational 



velocities in the disk agree. This procedure leads to an 
inclination angle i — 60° ± 5°. 

The blueshifted system B is most likely explained by 
the presence of outflowing gas. At various p. a. blueshifted 
structure in the wings of [Om] were earlier described by 
Davies & Illingworth (1986). Although they had not at- 
tempted a quantitative numerical decomposition of their 
line profiles, their insightful discussion led to the detection 
of widespread outflow on the eastern side of the nucleus. 

The spatial covering of the present data is insuffi- 
cient for clues about specific outflow geometries of com- 
ponent B. Also, prior to deeper interpretation the degree 
of uniqueness of the fit-procedure has to be discussed. To 
this end, we had also attempted fits with other constraints, 
e.g. fixing the blueshift of B (at 400 km/s relative to .4) 
rather than fixing the line width. The kinematics of the ro- 
tational component A was hardly altered suggesting that 
common velocities and linewidths in lines from different 
ions are a robust feature for a bulk component. 

Despite general agreement, to the northeast Table |] 
shows differences in the outflow velocities of up to ~ 
70 km/s. Such differences are within the accuracy of the 



fitting procedure that employs a simple basic function. 
Note that most differences are significantly below the com- 
ponent widths and the spectral resolution. The worst dis- 
crepancies are at 0'.'68 southwest, but here the outflow 
component starts to vanish. 

Row 'mean- 1502'inTable| gives the l.o.s. components 
of the outflow velocity relative to the emissivity maxi- 
mum, identified at 1502 km/s. On the northeastern side it 
is blueshifted by ~ (450 ±70) km/s; at 0768 southwest less 
consistent fits yield a blueshift of ~ 200 km/s. 

With a disk inclination angle of 60° the northeast out- 
flow velocity would be ~ 900 km/s relative to the emis- 
sivity maximum if it were directed along the minor axis. 
However, a detailed deprojection has to await more com- 
plete data. 



4.4. The Seyfert Mrkl210 
4.4.1. The WR-feature 

The three spectra extracted from F116 (Table 1) with 
2" wide bins from the central 6" (diameter) exhibit 
HeIIA4686 with an observed intensity of 25% of H/?. 
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Concentrated in the central 2" we detect a Wolf-Rayct 
feature around A4686 (for Mrk 1210 first shown in Storchi- 
Bergmann et al. 1998). 

Considering individual features, we note that 
[Fem]A4658.10 is clearly present, [Ariv]A4740.34 is 
present outside the innermost 2", while [Anv]A4711.34 
and [Neiv] A4724.15 are possibly present, but cannot be 
clearly discerned at our signal-to-noise level. 

The WR-feature is an unambiguous sign for the pres- 
ence of a starburst which appears to arise cospatial with 
the NLR or in a region roughly seen in projection to the 
inner kiloparsec in which the bright part of the AGN NLR 
is situated. 




4.4.2. Analysis based on single narrow-line component 

Firstly, peak positions of the strongest lines were mea- 
sured. The derived nuclear velocity curve along p. a. 
90° yields a small velocity gradient (Fig. [l5]). The gradi- 
ent is flattened because of seeing and guiding inaccuracies 
of ~ 2". A deconvolution is not warranted because of the 
lack of an appropriate point-spread function. 

Locating the kinematical center KC midway within the 
slope yields a measured difference to the brightness maxi- 
mum BM of 1'.'5. In the weaker blue spectrum (F116) the 
slope plus three extraction rows westwards can be mea- 
sured in [Oiii]. Here the difference between KC and the 
[Om] BM amounts to two rows or 1'.'3. Due to the spatial 
smearing the real difference is likely ~ 1" or slightly less. 

From velocities measured in the weak spectra (F085, 
F087), which were obtained offset from the brightness cen- 
ter by 2" (slit along p. a. 90°) we estimate a p. a. of the 
kinematical major axis of 110° ± 20°. The uncertainty of 
the p. a. is due to the faintness of the extranuclear lines. 

From the observed slight velocity variations one cannot 
safely distinguish between bipolar outflow and rotation. 
We tend to regard it as a reflex from the rotating nuclear 
kpc-spiral disk seen in a snapshot-HST-image (Malkan et 
al. 1998). However, neither from the kinematical data nor 
from the printed HST image a meaningful value of incli- 
nation can be derived, so that no reliable mass estimate 
can be given. 

Leaving aside the global kinematics, various trial fits 
of the strongest line profiles (observed to the west of the 
location of the velocity gradient) were undertaken which 
led to reasonably successful single-component fits. For jus- 
tification of the results given in Table^ and, in particular 
for an alternative decomposition of the profiles described 
in the next subsection, some details of the procedures em- 
ployed are to be explained. 

Firstly, it was attempted to fit the apparent three-line 
blends Ha+[Nn] by either a sum of three Gaussians or 
three Lorentzian functions only constrained by the 3-to-l 
intensity ratio (and same width) of [Nil] . These fits badly 
failed, although the Lorentzians yielded smaller but yet 
unacceptably large residuals. The ratio of the measured 
peaks of A6583 and A6548 in the blend is only 1.8 to 1 



Fig. 15. Velocity curve (heliocentric) of Mrk 1210 based 
on the mean of the peak positions of Ha and [Nn]A6583 
in the raw spectrum F066 along p. a. 90° (east to the left). 
The brightness maximum is at the zero position. The curve 
is smoothed by spatial smearing (> 2"). 



instead of the expected 3-to-l ratio. This indicates that 
the flux from an additional component is stronger beneath 
A6548 than below A6583, which can be relatively easily 
achieved by the contribution of a broad Ha component in 
addition to narrow Ha. The wide Lorentz wings of a single 
Ha component (allowed to be broader than the [Nil] lines) 
cannot lift A6548 sufficiently. 

All trial fits with an additional broad component as- 
sume it to be at the same central wavelength as the nar- 
row Ha feature. One broad component, either Gaussian 
or Lorentzian, does not lead to a satisfactory fit if three 
narrow Gaussians are added, due to their typical 'lack of 
wings'. Solely equipped with Gaussians, one would need 
many more than these four components for an acceptable 
fit. 

Three narrow and one broad Lorentzian function pro- 
duce too much flux in the outskirts of the blend suggest- 
ing that the real underlying broad component does not 
have extended wings. Indeed, the best four-line fits were 
obtained with three narrow Lorentzians and one broad 
Gaussian concentric with the narrow Lorentzian of Ha 
(as shown in Fig. p~6|) . 

Qualitatively, the presence of a broad Ha component 
is not unexpected in view of the polarized (degree of polar- 
ization ~ 1%) broad Ha found in Mrk 1210 (Tran 1995). 
However, the faint polarized flux should not give such a 
strong signal in total flux. 

The good single-component fits of [Nil] suggest to fit 
the other strong forbidden-systems by single Lorentzian 
functions as well (with the joint constraints as given in 
Appendix |b|). The resulting velocities and line widths from 
these single-component fits are given in Table^. In general, 
these simple fits reproduce the gross structures relatively 
well and may be considered as a success for the Lorentzian 
functions because single Gaussians fail to represent the 
wings of all major lines. 
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Table 6. Heliocentric velocities and line widths (FWHM) (separated by | ) in km/s along p. a. 90° through the nucleus 
of Mrkl210 as derived from single-component Lorentzian fits. The zero position is that of the interpolated brightness 
maximum (symmetry center of brightness). 



Position 


H/3 


[Om] 


[Oi] 


Habroad 


iiLt narrow 


[Nn] 


[Sil] 


3'.' 02 E 






3992|809 


3999|2268 : 


3999|116 


3998| 137 


3999| 148 


2'.' 34 E 


4008| 329 


4001 1471 


4076 715 


4004 1893 


4004 188 


4019213 


4017 177 


1766 E 


4021 450 


4034 507 


4104 605 


4048 1759 


4048|313 


4052 296 


4053 256 


0798 E 


4056|557 


4054|495 


4127|617 


4076| 1733 


4076| 331 


4081|287 


4072 1 274 


0730 E 


4066|595 


4063 1 494 


4133 607 


4088 1707 


4088|351 


4092 319 


4081 276 


0738 W 


4043 486 


4051 572 


4135 587 


4089 1703 


4089| 392 


4095 328 


4083 271 


1706 W 


4055 474 


4058|505 


4126 631 


4084 1684 


4084|396 


4091 319 


4084275 


1774 W 


4017|421 


4042 1 478 


4097|640 


4078|2429 


4078| 267 


4074|212 


4088|227 


2742 W 








4092 5129 : 


4092|137 


4076 144 


4094 166 




6550 6600 6650 6700 



VAngstr. 

Fig. 16. Ha+[Nn]AA6548,6583 from 0768 cast of the re- 
gion where lines are at maximum. The total fit (dashed 
line) consisting of a sum of three Lorentzians (continu- 
ous line) for narrow Ha and [Nil] (the [Nil] lines are con- 
strained by the 3:l-ratio and the same width) plus a broad 
Gaussian (dotted) concentric with the narrow Lorentzian 
that fits the core of Ha. 



Even H/3 can be represented by a single Lorentzian, 
somewhat wider than the above fit of narrow Ha and much 
narrower than broad Ha (Fig.^, Table g). This means 
that with the chosen fitting procedure the profiles of Ha 
and H/3 do not agree. The broad+narrow component de- 
composition of Ha cannot be transferred to H/3 by allow- 
ing for different values of extinction for the components. 

A close look unveils marked imperfections in the 
forbidden- line fits as well. E.g., the [Om] lines are asym- 
metric to the blue (Fig.|l7|), in the [Sn] blend a resid- 
ual excess remains in the red wing (Fig.[l8]). [Oi] shows 
structure and is redshiftcd relative to the other lines 
(Fig.||, Table§). In addition, the [Om] and [Oi] lines 
(w ~ 600 km/s) are significantly wider than the [Nil] and 
[Sn] lines (w ~ 300 km/s). Altogether, a more sophisti- 
cated decomposition is required which will be proposed in 
the next subsection. 
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Fig. 17. As in Fig.|T|, but for H/3 and [Om] AA4959, 5007 
as fit by single Lorentzian functions (dashed). The [Om] 
lines are constrained by the 3:l-ratio and equal width. 
Obviously they have a blueward wing asymmetry. The 
smooth continuous line shows the analogous, but much 
inferior fit with Gaussians. 

4.4.3. The decomposition of the line profiles 

Table || shows that the velocities of the single-component 
fits for emission from between 07 98 east and 1706 west 
agree within 10 km/s so that we combined the lines to 
represent the total emission from the ±1" of the bright- 
ness maximum. This gives practically no difference to fits 
of single extraction-rows, but reduces the noise and is 
more appropriate in view of the width of the spatial point- 
spread function of ~ 2". 

The red wings seen in the [Sil] and [Oi] profiles and 
the large widths of [Oi], H/3 and [Om] suggest the pres- 
ence of different kinematical components. However, the 
H/3 and [Om] profiles and the Ha+[Nn] blend are lacking 
clear substructure for providing anchoring features to con- 
strain different components. At first, we therefore looked 
for guidance by visible structures in other lines and at- 
tempted to apply the resulting components with the same 
widths and positions to the smoother profiles. 
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Fig. 18. [Sn]AA6716, 6731 from the same region as in 
Fig.|l6| as fit by two Lorentzians (dashed) with same 
width and position in velocity space. A fit with two 
Gaussians and otherwise same presumptions is displayed 
by the smooth continuous line. Athough not perfect, single 
Lorentzians appear to be better suited as basic fit func- 
tions than single Gaussians. 



2 - Mrk 1210, [01] 



^V^Lr-.;v^;~- 
6300 




Fig. 19. As in Fig.|§ but for [Oi]AA6300, 6363 as fit by 
two Lorentzians (dashed) with same width and position 
in velocity space. Smooth continuous line: The analogous 
trial fit with Gaussians. 



Fig. |o| shows a three-component Lorentzian fit to [Sn] . 
M. (not displayed) is the 'main' component determined by 
the peaks of the profiles, which therefore has the closest 
relation to the velocity curve in Fig.|l5|. The red-shifted 
component 1Z has been invoked to fit a residual shoul- 
der on the red side and the inter-line regime. The relative 
strength of the two [Sn] lines of A 7 ! is a free parameter. As 
in the case of NGC 1052, we assume that the particular 
strength of [Oi] is partially due to the second component 
(here: TV) with n 3> 10 3 cm~ 3 , which fixes its intensity 
ratio [Sii]A6731/[Sii]A6716 = 2.3. The components are M 
('main') and 1Z ('red') defined by (v=heliocentric veloc- 
ity, w=FWHM), both measured in km/s. The fit gives 
(v,w)(M) = (4067, 208) and (v,w)(TZ) = (4304, 350). 



Table 7. Relative intensity ratios of components Ai, 1Z 
and B as defined by corresponding velocity-width param- 
eters (v,w)(M) = (4064, 208), (v,w)(lZ) = (4304, 350) 
and (v,w)(B) = (3895, 727) (in km/s). The intensity of 
the line-core component M is set to unity. (1) and (2) cor- 
respond to alternative fits of the Ha+[Nll] blend (Figs. [2^ 
and GB) , respectively. 



Line 


K(l) 


6(1) 


U{2) 


B(2) 


[Oi] 


1.65 


2.05 






[Sn]A6716 


0.18 


0.07 






[Sn]A6731 


0.41 


0.17 






[Nil] 


0.42 


1.74 


0.29 


0.90 


Ha 


0.42 


1.73 


0.76 


1.71 


H/3 


0.65 


2.24 






[Qui] 


0.44 


1.97 







The weak third component B slightly improves the fit of 
[Sll], but its presence is inferred from [Oi] and [Om]. 

Applying M. and 7Z to [Oi] and [Om] takes care of 
the shoulder in the upper red wing of [Oi] and yields an 
excellent fit of the red wing of [Om], but, for the complete 
profiles, an additional 'blue' component B = (3895, 727) 
is also required. B is an excellent supplement to fit both 

23) 



the blue parts of [Oi] an [Om] (Figs.|2l]and |22|). Only the 
amplitudes of the components are left as free parameters 
to be determined during the fitting procedure; widths and 
positions are the same in all lines. Also H/3, ten times 
fainter than [Om]A5007 and therefore noisier, could be 
easily fit by a sum of M, B and TZ (Fig. ^3]). 

To find a corresponding fit of Ha+[Nn] with six free 
parameters (the amplitudes of Ha and [Nil] for each sys- 
tem) turned out to be too demanding for the software. 
Setting manual constraints while guided by the amplitude 
ratios of the other fits and the response to parameter vari- 
ations led to two acceptable fits, dubbed solution (1) and 
(2) (Figs.|4[and |§ Tables | and |) . 

In Fig. p3 the main-component amplitudes h(A4) of 
Ha and [Nn]A6548 are left free, while, in a series of fits, 
the corresponding h(7Z) and h(B) were manually adjusted, 
which arc, however, constrained to have fixed ratios to 
h(M). With this procedure each component is confined 
to the same [Nil] /Ha ratio. We obtain [Nll]A6583/Ha = 
0.538 and h{M):h{1Z):h{B) = 1 : 0.25 : 0.50 (solution (1)). 

However, to fix [Nil] A6583/Ha for all components is a 
much too simplifying assumption since it depends on den- 
sity and input flux in a more complicated manner than, 
e.g., the Balmer normalized intensities of [Oi] or [Om] (e.g. 
Komossa & Schulz 1997). [Nil] /Ha might differ from com- 
ponent to component. Therefore, in another series of fits 
different values of the relative [Nil] amplitudes h(JZ):h{M) 
and h{B):h{M) were predefined while the Ha amplitudes 
of the components were left as free fit parameters. The 
best fit (solution (2)) is shown in Fig. gq. The resulting 
intensity ratios are given in Tables and |S[ 
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Due to a lower x 2 we have a formal reason to pre- 
fer solution (2). The second, more convincing, reason is 
the more plausible physical assumption underlying fit (2). 
To achieve uniqueness in fitting the Ha+[Nll] blend more 
constraints than just fitting profiles have to be taken into 
account. 

In solution (2) [Nn]A6583 is slightly above the observed 
level. The discrepancy cannot be removed with redward- 
asymmetric components. In this regard, a broad compo- 
nent would be more flexible (see Fig.[l6|). Given the clear 
evidence for broad Ha in polarized flux (Tran 1995), a 
faint broad component underlies the lines in total flux as 
well. The fits shown here nevertheless demonstrate that it 
is not necessary to postulate a broad component. 

Table || gives the major line-intensity ratios for 
the adopted decomposition, which all lie in the AGN 
regime of the standard diagnostic diagrams of VO. The 
[Om]A5007/H^ ratios lie in the range 10 ± 2. Tl is char- 
acterized by relatively strong [Oi]/Ha combined with rel- 
atively weak [Nil] /Ha and [Sn]/Ha. 

We carried out further experiments to check the 
uniqueness of the fits. It turned out, that the visible nar- 
row cores of the lines constrain A4 rather well, and next 
comes 1Z determined to fill in the shoulders in [Sn] and 
[Oi]. The weakest constraint is for B, which takes account 
'for the rest' and could be widened to ~ 900km/s (and 
correspondingly shifted) on the expense of the other com- 
ponents. However, experiments in this direction yielded 
worse total fits so that for the final discussion solution (2) 
given above is preferred. 
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Fig. 20. [Sll]AA6716, 6731 from the central 2" of Mrk 1210 
with maximum line emission as fit by three kincmatically 
separated Lorentzian emission systems: a blue component 
B (dotted), which is unusually faint in [Sil] ; a red compo- 
nent 1Z (dash-dotted) (explaining the redward asymmetry 
of [Sn] ; the main component M. that fills the cores of the 
lines (not shown)). The dashed line shows the total fit 
made up of M + B + TZ. 
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Fig. 21. As in Fig.[|| but for [Oi] AA6300, 6363 from the 
maximal line-emitting 2" of Mrk 1210. 



12 


- 




e 


Mrk 1210, [OIII] 1 


1 


4 




J 
ft 














5000 5050 5100 

\/Angstr. 



Fig. 22. As in Fig.|o|, but for [Om]AA4959, 5007 from the 
maximal line-emitting 2" of Mrk 1210. 

5. Discussion 

5.1. The narrow-line profiles of NGC 1052 and 
Mrk 1210 

5.1.1. Kinematic bulk components 

The slightly novel ansatz of the present work lies in the 
attempt to disentangle two to three emission-line systems 
from the observed line profiles. Ab initio such a task may 
be subject to pitfalls like non-uniqueness and particular 
care has to be devoted to the choice of a basic function 
that represents a single profile. 

The above fitting experiments carried out with the 
spectra of NGC 1052 and Mrk 1210 showed that single- 
component fits as well as multi-component fits led to lower 
residuals and a lower number of fitting components when 
using Lorentzians as basic functions rather than the com- 
monly employed Gaussians. The additional surprising ef- 
fect of this choice was that we could essentially dispense 
with broad Ha (and H/3) components. There is no doubt 
that both galaxies display broad Ha components in polar- 
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Table 8. Logarithmic intensity ratios as used in the diagnostic diagrams by Veilleux & Osterbrock (1987; VO) of 
components M, TZ and B of Mrkl210. (1) and (2) correspond to alternative fits of the Ha+[Nll] blend (Figs. ^4] and 



Line 


M(l) 


11(1) 


B(l) 


M(2) 


11(2) 


B(2) 


log([Om]A5007/H/3) 


1.075 


0.906 


1.019 


1.075 


0.906 


1.019 


log([Nn]A6583/Ha) 


-0.269 


-0.269 


-0.270 


-0.107 


-0.529 


-0.385 


log([Sn]A6716/Ha) 


-0.515 


-0.879 


-1.883 


-0.490 


-1.109 


-1.854 


log([Sn]A6731/Ha) 


-0.504 


-0.518 


-1.523 


-0.480 


-0.747 


-1.491 


log([Sn]AA6716,6731/Ha) 


-0.211 


-0.361 


-1.366 


-0.184 


-0.591 


-1.334 


log([0 i]A6300/Ha) 


-0.647 


-0.055 


-0.574 


-0.343 


-0.285 


-0.543 
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Fig. 23. As in Fig.|o|, but for H/3 from the maximal line- 
emitting 2" of Mrk 1210. 
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Fig. 24. As in Fig.[2(], but for Ha+[Nll] from the maxi- 
mal line-emitting 2" of Mrk 1210. In this special fit, the 
[Nil] /Ha ratio is the same for systems M, B and TZ (so- 
lution (1), see text for more details). 

ized flux but these are faint and are unlikely to be visible 
in spectra taken without a polarizing device. 

If the Lorentzian decomposition is correct, the broad 
Ha component extracted by Ho et al. (1997) in the total- 
flux spectrum of NGC 1052 has to be considered as spuri- 
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Fig. 25. As in Fig.|24|, but here the intensity ratios 
[Nil] /Ha are allowed to differ among systems A4, B and 
TZ (solution (2), see text for more details). 

ous. Our alternative decomposition at least suggests to be 
cautious with claiming a broad Ha component hidden in 
the Ha+[Nn] blend if it is only based on a mere decompo- 
sition in total flux. In this view it seems to be a coincidence 
that spectropolarimetry finally confirmed NGC 1052 as a 
LINER- 1. 

Before giving possible astrophysical settings for 
Lorentzian profiles, an empirical justification is in or- 
der. The significance of the decomposition into different 
Lorentzian components rests on the internal consistency 
achieved for all strong line profiles from H/3 to [Sn]A6731. 
In particular, the same decomposition was employed for 
the collisionally excited transitions of [Oi] , [Oni] , [Nil] and 
[Sll] as for the allowed recombination lines Ha and H/3. 
Because of electron densities believed to exceed 10 9 cm~ 3 , 
a classical AGN BLR would only exhibit broad allowed 
lines (with line widths (FWHM) in the range 10 3 ~ 4 km/s) 
so that we consider the components extracted here (with 
FWHM ~ 10 2_3 km/s) as belonging to parts of a NLR 
distinguished by different bulk motions. 

Density tracers usually suggest densities in the range 
10 2 ~ 4 cm~ 3 for the main components of NLRs. For 
NGC 1052 and Mrk 1210 blue components B with indi- 
cations for rather high densities (~ 10 6 cm" 3 or more) 
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were isolated. In particular, B of Mrkl210 hardly shows 
up in [Sn]AA6716, 6731, which suggests n e ^> 10 3 cm~ 3 . 
However, since the critical densities of the well observed 
doublet [Nn]AA6548, 6583 are ~ 10 5 cm" 3 , the latter value 
should not be exceeded by several orders of magnitude 
Furthermore, the trends shown in Komossa & Schulz 
(1997) indicate, that to boost [Nil] /Ha relative to [Sii]/Ha 
at high densities, a large ionizing input flux is necessary. 
B might therefore be an 'ILR' (intermediate-line region) 
with a density and a distance intermediate between the 
BLR and NLR. This designation is also supported by its 
large line width of 727km/s, somewhat intermediate be- 
tween a 'typical' NLR and a typical BLR. 

In the early days of NLR analysis width-variations of 
narrow lines were frequently attributed to stratifications 
in velocity space observationally manifested via different 
critical densities or ionization potential (e.g. Pelat et al. 
1981, Filippenko & Halpern 1984, Whittle 1985b, Schulz 
1987). The decompositions proposed here explain different 
line widths by summing line components, each of constant 
width in velocity space, but different weighting for a com- 
ponent in each line. This assumption implies that each 
emission-line cloud emits all lines considered here, which 
will be the case when they are mainly ionization bounded. 

While our data suggest that the M. components arise 
in a rotating disk with 'normal' densities (their velocities 
fit the rotation curves), the TZ and B components appear 
to arise in bulk outflow systems. 



FWHM = 20. 
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Fig. 26. A Lorentzian (dashed) and a Gaussian (dotted) 
with same full width at half maximum. 



5.1.2. Why Lorentzians? 

Gaussians are relatively well suited to describe instru- 
mental functions or thermal Doppler broadening or some 
kind of turbulence induced by random motions with 
Maxwellian-type velocity distributions. Lorentzians differ 
from Gaussians with same FWHM by a marginally nar- 
rower core and a slower fall-off to large velocities leading 
to more extended line wings (Fig.p6h. The non-Gaussian 



nature of observed [Oin] line profiles was early described 
by Whittle (1985a), who had noted that (i) "almost all 
the sample have [Om] profiles with a stronger base rel- 
ative to the core than Gaussian" (his Sect. 5.3) and (ii) 
"the extreme wing widths (FWZI) usually lie in the range 
1000-3000 km s^ 1 " (his Sect. 5.2). 

For our line profiles, the observed line wings would 
have typically reached the noise levels at ~ 6% of the 
maximum. Eq. |B.l| implies that a Gaussian falls to 6.3% 
of the maximum at two HWHMs (half width at half max- 
imum) from the center while a Lorentzian has decreased 
to 5.9% at four HWHMs. For components M, TZ and 
B of Mrkl210 described in Table| we get 4 HWHMs 
= 2w = 416, 700, 1454km/s, respectively. How can such 
large velocities of ~ 10 2 ~ 3 km/s in the line wings be at- 
tained? 

Turbulence of this order within the emission-line clouds 
would be strongly supersonic in gas with an observed tem- 
perature of T ~ 10 4 K, in which the characteristic ther- 
mal velocity and the velocity of sound are of the order of 
~ lOkm/s. Hydrodynamic shocks would result in dissipa- 
tion of the gained energy, most likely by radiation, perhaps 
contributing to, but not dominating the ionization equi- 
librium that is essentially governed by radiative input. To 
uphold the velocities, a strong steady mechanical energy 
input would be required. 

Outflow velocities of the order of 10 3 km/s have fre- 
quently been seen in AGN absorption lines. Elvis (2000) 
reviewed the observations and proposed an AGN-type-I 
unification scheme. Cloudlets comoving in a bipolar wind 
system might lead to such line wings whose width would 
strongly depend on the viewing angle of the system. 

A less explored, albeit not far-fetched, mechanism to 
produce the line wings might be magneto-hydrodynamic 
waves that could mimick supersonic turbulence (Shu 
1992). 

A different type of alternative would be scattering of 
photons at fast particles. Dust particles are far too slow, 
but electrons in the cited photoionized clouds have ve- 
locities of 10 3_4 km/s. Scattering of a line at these elec- 
trons could cause broad line wings via redistribution in 
frequency space. In the early times of AGN spectroscopy, 
this type of scattering was proposed to explain the widths 
of lines from the BLR (Weymann 1970; Mathis 1970). In 
a hotter intercloud medium with T ~ 10 6 ~ 8 K, which 
may confine the emission-line clouds by its thermal pres- 
sure, electrons can attain velocities of 10 4 ~ 5 km/s. Shields 
& McKee (1981) suggested electron scattering in such a 
medium as an explanation for very broad Ha wings ob- 
served in the quasar B 340. Their case could be transferred 
to the Lorentzian wings, if one sees only the 'tip of the ice- 
berg' of hot-ISM scattering because otherwise the wings 
would have to be much wider. 

The fraction of scattered flux is roughly given 
by the electron-scattering optical depth r e and the 
covering fraction of the scattering medium. Effective 
column-densities of fully ionized gas of 2 10 22 cm" 2 
would be required to obtain electron-scattering opti- 
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cal depths on the percent level. Attributing such a 
column to a kpc-extended sphere of a hot intercloud 
medium would lead to a total bremsstrahlung lumi- 
nosity (Eq. 5.15b in Rybicki & Lightman 1979) of 
2.6 1 42 K/10cm" 3 )V( r / 1 ° 6K )( i ?/500pc) 3 erg/s. This 
is a lower limit because the estimate neglects line emission. 
Thermal soft-X-ray emission is estimated to be ~ 10 40 
erg/s for NGC 1052 (Weaver et al. 1999), and is measured 
to be - 10 41 erg/s in Mrkl210 (Awaki et al. 2000). We 
conclude that we do not see a sufficient amount of hot 
coronal gas to produce the required wings. However, a vi- 
able alternative would be scattering at the cooler (T ~ 10 4 
K) emission-line clouds with electron densities ~ 10 2-3 . 
With a scale R ss 10 2 pc and a geometrical scattering 
factor / w 0.3 one easily obtains a ratio of scattered lumi- 
nosity to input luminosity L s / L- m = fr e of a few percent, 
which suffices to explain th extended wings of the narrow 
lines. 

5.1.3. Relationship to spectropolarimetry 

Spectropolarimetry of NGC 1052 (Barth et al. 1999) and 
of Mrk 1210 (e.g. Tran 1995) provides direct evidence for 
scattering of the continuum and broad lines. 

In Mrk 1210 the degree of polarization P in the forbid- 
den lines decreases significantly (to less than 1%; to ~ 0% 
in the cores of [Oiii] ) compared to that of the continuum 
indicating that the main body of the narrow-line profiles is 
not scattered (this is the situation in many 'hidden-Seyfert- 
1' galaxies, which are thereby backing the unified model). 
In NGC 1052 there is hardly a decrease in the forbidden 
lines, suggesting that a major part of the forbidden-line 
profiles is scattered as well, rendering this object as an 
unusual 'hidden-LINER-1' object. How can these obser- 
vations be reconciled with the line-profile decomposition 
proposed here? We discuss the galaxies in turn. 

The average angle of polarization in Mrk 1210 is ~ 30° 
(Tran 1995). This is 80° off from the kinematical major 
axis of the central kpc-disk (cf. Table ||), whose geometry 
is not favorable to give rise to polarized radiation because 
of its near-face-on view. A suggestive explanation is found 
by the observed shift of the line-continuum maximum by 
~ 1" to the west of the kinematical center (see Fig.|l5|, 
where the zero point of the abscissa is at the brightness 
maximum), which is consistent with a view that the true 
AGN is obscured and is situated at the kinematical center. 
Via some escape route AGN radiation reaches the bright- 
ness center, thereby giving rise to the polarized scattered 
broad-line wings and the polarized part of the continuum. 
In the AGN standard view, the high-excitation Seyfert- 
2 characteristics of the brightness center would be traced 
back to photoionizing AGN radiation reaching this region. 
Because the AGN continuum is only seen in scattered light 
and the stellar continua from the starburst plus the older 
stellar populations certainly dominate, it is not surprising 
that the precise amount of visible 'featureless continuum' 
is under debate (Schmitt et al. 1999). 



In NGC 1052, Barth et al. (1999) found a polarimet- 
ric position angle of Ha (+[Nn]) of 178°, which is neither 
aligned with - nor perpendicular to - the major axis of 
the gaseous disk (p. a. ~ 45° following Davies & Illingworth 
(1986)). However, this apparently accreted narrow-line gas 
does not necessarily reflect the orientation of the 'central 
machine'. If the kpc-radio structure detected by Wrobel 
(1984) can be identified with a jet perpendicular to the 
central accretion disk, the polarimetric p. a. would per- 
fectly fit with the idea that the polarized component is 
scattered light originating in a nucleus obscured from di- 
rect view. 

However, unlike in Mrk 1210 (Tran 1995; second panel 
in his Fig. 15), in NGC 1052 the degree of polarization 
P (Barth et al. (1999); their Fig. 1, middle and lower 
panel, and Fig. 2, lower panel) only hardly decreases in 
the [Nll]+Ha blend and shows upturns in [Oi] and [Sii] . 
Considerable noise and possible Galactic foreground po- 
larization affecting the spectrum lets Barth et al. focus on 
[Oi]A6300, whose value of P lies significantly above that 
of the continuum. In particular, they note that [Oi] is nar- 
rower in polarized light and that the [Oi]/[Sii] flux ratio 
is larger in Stokes flux than in total flux. 

These latter findings resemble those we found for 
component B in NGC 1052 (cf. Table||, last two rows). 
Relative to A, in B usually [Oi] is stronger and [Sn] is 
weaker. The dense narrow-line gas of B could be (par- 
tially) hidden and its radiation be scattered like the broad 
polarized wings of Ha. 

Alternatively, Barth et al. (1999) suggested transmis- 
sion through a region of aligned dust grains as a possi- 
ble mechanism for the [Oi] polarization. Indeed, Davies 
& Illingworth (1986) noted a dusty zone on the eastern 
side of the nucleus where component B is predominantly 
detected. 

5.2. Relationship to H2O megamaser activity 
5.2.1. General aspects 

The A^1.3cm transition of H2O can be driven to mase by 
collisions, provided that the densities and temperatures 
are high enough (^510 7 cm -3 and JJ400K; see Sect. 1). In 
galactic sources of water maser emission, namely in star 
forming regions, the necessary conditions can be produced 
by shock waves which provide not only proper physical 
conditions (e.g. Elitzur 1992, 1995) but also greatly en- 
hanced H 2 abundances (e.g. Melnick et al. 2000). The 
unique association of H2O megamasers with AGN and 
the fact that not all of the megamaser emission can arise 
from unsaturated amplification of an intense nuclear ra- 
dio continuum (e.g. Greenhill et al. 1995a; AGN brightness 
temperatures can be much higher than Tb^SlO 4 K, encoun- 
tered in ultracompact galactic HII regions) indicates, how- 
ever, significant differences between galactic maser and 
extragalactic megamaser sources. Possible scenarios that 
may produce suitable physical environments for luminous 
megamaser emission from the so-called 'disk-masers' (see 
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Sect. 1) are (1) irradiation of dense warped disks by X- 
rays from the AGN (Neufeld et al. 1994), (2) shocks in 
spiral density waves of a self-gravitating disk (Maoz & 
McKee 1998) and (3) viscous dissipation within the accre- 
tion disk (Desch et al. 1998). For the so-called 'jet-masers', 
interaction of dense neutral gas with a radio jet is a viable 
scenario that may also yield temperatures of several 100 K 
and large H 2 abundances (Peck et al. 2001). 

5.2.2. Outflow as H2O maser trigger: winds of an AGN 
or a starburst? 

In view of the overall angular scales, i.e. a few milliarcsec- 
onds for intcrferomctric maps of 1.3 cm megamaser emis- 
sion and a few arcseconds for our study, a connection be- 
tween optical and radio data is difficult to establish. In 
elliptical galaxies the large scale dust lanes and nuclear 
jet orientation is correlated (e.g. Kotanyi & Ekers 1979; 
Mollcnhoff et al. 1992). In Seyfert galaxies, however, large 
scale disks and nuclear disks tend to be inclined with re- 
spect to each other (e.g. Ulvestad & Wilson 1984). In 
some galaxies, the nuclear jet may be oriented not too 
far from the plane of the parent galaxy and may thus in- 
teract with molecular clouds, forming jet-masers off the 
very nuclear region (see e.g. Gallimore et al. 2001 for the 
case of NGC1068). The near-to- face-on large scale view 
onto Mrkl210 (Malkan et al. 1998) is thus not arguing 
against high shielding column densities along the line-of- 
sight toward the nucleus. 

Radio jets may be a strong agent to drive maser action 
when hitting high-density clumps or entraining molecu- 
lar gas (see Elitzur 1995 for physical processes involved). 
However, optically detectable ionized outflowing material, 
e.g. a magnetized wind from an accretion disk (cf. Krolik 
1999) that would be less violent close to the equatorial 
plane than towards higher latitudes, may also be relevant. 
In systems with a nuclear thin disk and a thicker outer 
torus or in systems with misaligned angular momenta 
of nuclear and outer disk, a nuclear wind impinging on 
molecular clouds farther out appears to be an attractive 
alternative to the previously outlined jet-maser excitation 
scenario. As in the case of an interaction with a jet, shocks 
at speeds not destroying the dust grains and a warm (sev- 
eral 100 K) H2O enriched postshock medium are expected. 
The spatial extent of the outflowing gas that may have a 
larger 'cross section' than the faster radio jet makes such 
an interaction particularly likely. 

If one relies on outflow as one of the primary maser 
agents, it must be intriguing that H 2 megamasers have, 
until recently, only been detected in galaxies with AGN- 
typical features because large-scale winds and outflows are 
properties not only associated with AGN but also with 
starbursts. The key to this conundrum may be the fact 
that the sample of known H 2 megamaser sources is still 
rather small. Recently H 2 megamasers have been ob- 
served, for the first time, in starburst galaxies (Hagiwara 
et al. 2002; Peck et al. 2003) and for the 'traditional' ab- 



sence of H 2 megamasers in type-1 AGN, proposed by 
Braatz et al. (1997), a first counterexample may have been 
found as well (see Nagar et al. 2002). 

Below, some characteristic features of our four individ- 
ual galaxies are briefly summarized. 

5.2.3. Individual sources 

In NGC 1052, a notable feature is the presence of a 
strong high-density outflow (Sect. 4.3.2). For this galaxy 
the (VLBA+VLA) observed maser sources are found to 
coincide with the innermost part of the southwestern ra- 
dio jet rather than being related to the nearby putative 
nuclear disk (Claussen et al. 1998). The masers are red- 
shifted by 80 to 200km/s relative to v sys as tabulated by 
de Vaucouleurs et al. (1991). The distance to the 43 GHz 
radio core, the supposed location of the central engine, is 
less than a milliarcsecond (0.07 pc). 

Due to the spatial coincidence Claussen et al. (1998) 
suggested to relate the maser activity to the VLBI jet 
or to a fortuitously positioned foreground cloud that am- 
plifies the radio continuum of the inner jet. The strong 
high-density outflow disentangled from the rotating disk 
(Sect. 4. 3. 2) raises another alternative source to excite the 
maser activity. In this case the small solid angle of the 
inner jet, the presumably larger solid angle of the out- 
flowing component and the location of the masers in front 
of the brightest parts of the jet suggest unsaturated H 2 
emission amplifying the strong radio background. 

Towards Mrk 1210, the brightness maximum BM of the 
emission-line region is rather compact (Falcke et al. 1998) 
and BM is slightly misplaced from the kinematical center. 
Strong radial flows in BM may be related to the flows that 
give rise to megamaser activity. No interferometric radio 
map of the H 2 emission has yet been published. 

In IC2560 and NGC1386 we found little obvious evi- 
dence for outflow, presumably because of the near-edge-on 
view on these galaxies making a flow perpendicular to a 
disk hard to detect. The NGC 1386 line profiles show pro- 
nounced structure, and Rossa et al. (2000) isolated even 
more subtle substructures suggestive of non-circular mo- 
tions like locally expanding gas systems. Such more local- 
ized flows could be related to megamaser activity as well. 
However, we caution that integrating over l.o.s., which in- 
tersect a dusty and patchy rotating disk, can well lead 
to bumpy line profiles. Interpreting such bumps by lo- 
cal Gaussians may be misleading. A meaningful discussion 
requires observations with significantly higher resolution 
and a detailed model of the dusty spiral disk. While the 
megamaser in IC 2560 may originate from a circumnuclear 
disk not related to outflowing gas (Ishihara et al. 2001), 
no high resolution H 2 map has yet been published from 
NGC 1386. 

In passing, it should be noted that the mere pres- 
ence of outflow in H 2 megamaser galaxies is not new. In 
much more detail outflow bulk systems were studied in the 
nearby galaxies M 51 (Cecil 1988) and NGC 4258 (Cecil et 
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al. 2000), but a direct relationship to their megamaser ac- 
tivity (e.g. Herrnstein et al. 1999; Hagiwara et al. 2001) 
has so far been elusive. M51 and NGC4258 are likely not 
suitable candidates to prove such a relationship, but op- 
tical data with subarcsecond resolution are needed for a 
convincing demonstration in other sources. 

6. Final conclusions 

We have analyzed exploratory optical spectra of four H2O 
megamaser galaxies in order to disentangle kinematically 
discernable bulk components of their emission- line regions. 
For each object we propose a model of the global kine- 
matics. The gas dynamics in IC 2560 and NGC 1386 is 
mainly determined by gravitationally dominated motions 
in a disk. However, there is a bias due to the close-to- 
edge-on view onto their disks making minor-axis outflow 
hard to detect. For NGC 1052 and Mrk 1210 emission from 
gravity-dominated dynamics of orbiting clouds could be 
separated from that of systems being in outflow. Line- 
broadening within the individual bulk components ap- 
pears to be characterized by scattering wings. 

IC 2560 is a highly inclined Sb spiral with a nuclear 
high-ionization Seyfert-2 spectrum displaying relatively 
narrow (FWHM ~ 200km/s on an arcsecond scale) emis- 
sion lines from a central rotating disk. If its major-axis 
orientation and inclination agree with that of the outer 
galaxy, the core mass (radius 100 pc) will be 91O 6 M0, 
consistently larger than the 2.81O 6 M , given by the or- 
bital velocities of the maser sources presumably embedded 
in a compact disk of outer radius 0.26 pc as measured by 
Ishihara et al. (2001). A slight blueward asymmetry of the 
line profiles suggests the presence of outflow (Sect. 4.4.1, 
Fig. 3). 

The nuclear mass of - 1O 7 M in IC2560 appears to 
be rather small. It would be larger by a factor of 1.5 if a 
distance of 39 Mpc (see comment in Table 2) were adopted. 

NGC 1386 is a highly inclined Sa spiral with a nu- 
clear high-ionization Seyfert-2 spectrum displaying rather 
structured asymmetric emission-line profiles that we inter- 
pret as arising from various parts of a central edge-on ori- 
entated disk. We dismiss the term 'radiation cone' for this 
object because the elongated Ha morphology may simply 
reflect the edge-on view on the matter-bounded disk. The 
morphology together with analogy reasoning suggests that 
the disk is a dusty gaseous mini-spiral, which appears to 
be warped. Its rotation curve inside a diameter of 1.6 kpc 
yields a dynamically acting mass of ~ 5 10 9 M Q , while the 
visible line-emitting ionized gas has a mass of ~ 10 5 ~ 6 M , 
which is a very small fraction of the total mass of nuclear 
gas that is likely to be of the order of ten percent of the 
dynamical mass. 

There is a lack of line emission from the kinematical 
center (the center of the disk), which is probably due to 
obscuration by an equatorial dust disk that is evidenced 
by a dust lane crossing the center. A slight east- west gra- 
dient in velocity arising in protrusions seen on an HST 
image (Ferruit et al. 2000) could be due to minor-axis 



outflow, i.e. a flow essentially perpendicular to the disk 
(Sect. 4.2.2). 

AGN narrow-line emitting disks like in IC2560 and 
NGC 1386 explain general findings that the widths of 
narrow-lines are broadly correlated with the nuclear stel- 
lar velocity dispersion (Nelson & Whittle 1996). Such 
NLRs are likely to be part of the nuclear gas-dust spi- 
rals commonly seen on HST images of AGN on a scale 
of ~ 10 3 pc (Malkan et al. 1998; Martini & Pogge 1999; 
Martini 2001). However, so far unobserved very narrow 
lines should be emitted from galaxies with central disks 
seen face-on, if NLRs were only dominated by disk dy- 
namics. The other two (non-edge-on) galaxies provide ex- 
amples for additional dynamical components. 

NGC 1052 is an elliptical galaxy with a nuclear low- 
ionization emission line region (LINER). All major line 
profiles from the brightness center out to 2" northeast 
can be decomposed into two components A and B, each 
with the same width in velocity space for all lines at a 
given spatial position. A is attributed to a rotating disk, 
while B is blueshifted by ~ 400km/s suggesting outflow 
of high-density gas relative to the disk. Line ratios of both 
A and B are LINER-like. We propose to relate B to the 
peculiar forbidden-line polarized component detected by 
Barth et al. (1999) (Sect. pM.3| ). A faint broad Ha com- 
ponent, which is seen in polarized flux, is not required in 
total flux. The presence of strong high-density outflow sug- 
gests to explore the possibility whether such a flow could 
be an alternate trigger of the H2O megamaser activity in 
NGC 1052. 

Mrk 1210 is an amorphously looking Seyfert-2 galaxy 
that contains a nuclear kiloparsec-spiral seen face-on. A 
bright emission-line spot appears to be slightly shifted 
from the partially obscured kinematical center. In total 
(unpolarizcd) flux the emission from the bright spot can 
be decomposed into three narrow-line components M. (the 
'normal' main component), 1Z (a redshifted component), 
and B (a blueshifted component) without any BLR com- 
ponent. M. is tentatively attributed to the rotating spiral 
from which only a small velocity variation could be mea- 
sured because of its near-face-on orientation. TZ and B 
are outflow components that may both be related to the 
cospatial starburst and to ionization by the active nucleus. 
The megamaser activity might be related to the outflow. 

The line ratios from the bulk components of Mrk 1210 
fall into the AGN regime of the VO diagrams. The shifts 
of component 1Z (by 240km/s to the red of M) and 
of component B (by 169km/s to the blue relative to 
Ai) are moderate, but of significant influence because 
major portions of the emission-line region are involved. 
Regarding line widths and densities B shows characteris- 
tics of an intermediate- line-region between BLR and NLR 
(Sect. pUl| ). 

Finally, it should be noted that the successful line- 
profile decompositions rest on Lorentzians rather than 
Gaussians as basic functions. Hence, extended wings of in- 
trinsic narrow-line profiles of the bulk components may not 
be uncommon. Crude estimates suggest that electron scat- 
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tcring at the ionized gas itself might lead to a viable ex- 
planation, but detailed modelling of such processes would 
be useful. 

There exists a possible connection between optical and 
radio lines. It is suggested that nuclear outflows impinging 
onto dense molecular clouds may provide a suitable trig- 
ger for megamaser emission. If inner torus and outer disk 
are misaligned, such an interaction appears to be partic- 
ularly likely. A convincing verification of such a scenario 
requires, however, optical spectroscopy with subarcsecond 
resolution. 

Appendix A: Extraction of emission lines 

A line intensity is obtained by integrating over a line pro- 
file in a flux calibrated spectrum. Both this integration 
and any profile decomposition are usually carried out af- 
ter having subtracted the template spectrum and, if this 
has not led to a 'zero base line', a normalizing continuum 
has been subtracted as well. The latter is defined by a 
spline fit through marks set with the mouse cursor. The 
subjectivity of this eye fit hardly affects our basic results 
which rely on strong lines. The second step is necessary in 
case of an inappropriate template and/or the presence of 
a nonstellar continuum. 

For each object, we define the spatial 'center' of an 
emission line region by the intensity maximum of the emis- 
sion lines. In all cases this center coincides with the max- 
imum of the continuum between the lines within a small 
fraction of an arcsecond. This kind of center may not co- 
incide with the kinematical center (usually the 'nucleus' 
of a galaxy) derived from an assumed symmetry of the 
velocity curve. 

In the main text, discussion of emission- line intensity 
ratios is largely confined to the diagnostic classifiers pro- 
posed by VO which do not require large and notoriously 
uncertain reddening corrections. Occasionally intensity ra- 
tios of the [Sn] pair A6716 and A6731 will be given, which 
may be converted into some kind of weighted mean of the 
electron densities in the emission-line gas via Fig. 5.3 in 
Osterbrock (1989). 

Appendix B: Common line-profile analysis 
procedures 

In AGN even the 'narrow' lines show considerable profile 
structure if measured at sufficient resolution (Vrtilek & 
Carleton 1985; Whittle 1985a). When single-component 
Gaussian fits fail, commonly multi-component Gaussian 
fits arc employed. However, our multi-component fitting 
experiments of the spectra from NGC 1052 and Mrk 1210 
showed that in these cases Lorentzian functions C(AA) 
(also called 'Lorentzians'; C stands for 'Cauchy func- 
tion') appear to be better suited if one intends, apply- 
ing Occam's razor, to minimize the number of required 
components and parameters in a multi- function fit. This 
may not be generally true. An instructive case of an asym- 
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metric line profile being Lorentzian on the blue side and 
Gaussian on the red side is given by IC 2560 (cf. Sect. 4.1). 

In a crude picture, the basic kinematical configurations 
expected are rotating disks, cones of outflowing gas or ex- 
panding shells. As shown in Schulz et al. (1995), line-of- 
sight integrations through such configurations would nei- 
ther yield Gaussians nor Lorentzians. However, Schulz et 
al. (1995) also showed that the spatial and spectral smear- 
ing by an instrumental function that is wide in comparison 
to kinematical gradients and the presence of strong ran- 
dom motions superposed on the systematic bulk motions 
will lead to profiles characterized by a more symmetric 
core- wing structure. 

This case applies here because our spectral resolution 
of <~ lOOkm/s and spatial resolution of <~ (1 — 2)" leads 
to integration over relatively large cells of the phase space 
of a NLR. Hence, for finding bulk motion components ex- 
perimental fits are carried out by making use of Gaussians 
and Lorentzians as basic functions. 

Lorentzian and Gaussian functions are fixed in the 
MIDAS data reduction software by the three convenient 
parameters central line position Ao (so that A A = A — A ), 
w = FWHM = full width at half-maximum and h, the 
maximal value of the profile function: 

C(AA) = A \ 9/ „ (B.l) 

v ; l + 4AA 2 /w 2 v ; 

G(AA) = fcexp(-ln(2) (4A\ 2 /w 2 )) (B.2) 

The ratio of the intensities of lines arising from 
the same upper level is given by the ratio of 
the corresponding transition probabilities (taken from 
Tables 3.8 and 3.10 in Osterbrock (1989)). Therefore 
the fits of [Om]AA4959,5007, [Oi]AA6363, 6300 and 
[Nn]AA6548, 6583 are constrained by fixing intensity ra- 
tios in the proportion 1:3 and by equalizing the FWHM 
w. 

Aside from these constraints, the basic strategy of fit- 
ting the line profiles is to start with a minimal set of as- 
sumptions. At first it is attempted to fit each narrow line 
by a single Lorentzian function or a Gaussian. If this fails 
we look for features (shoulder, second peak, asymmetric 
wing) giving clues for the initial condition to start a fit 
with a second component. In case of success, this compo- 
nent (with same velocity width w, where (w in km/s) = (c 
in km/s) x (FWHM in A) / (Ai ine in A)) will be applied to 
the other strong lines. If the positions among different lines 
agree sufficiently well in velocity space, the fit is consid- 
ered as a success. In case of Mrk 1210 a three-component 
fit turns out to be necessary for which the positions of 
the components are also fixed only leaving the component 
intensity as a free parameter in different lines. 

In physical terms, a bulk-motion system would emit 
lines of approximately the same width if it is, e.g., com- 
posed of local emitters (with sharp lines compared to 
the spectral resolution), which are photoionized to such 
a depth that they emit all lines used in the analysis. 

Component line widths (FWHM) will be given as mea- 
sured, i.e. without correcting for our finite spectral reso- 
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lution of lOOkm/s. Applying the convolution theorem for 
two Gaussians, a measured FWHM of 200km/s would be 
14% too large, with a rapidly decreasing relative excess 
for larger FWHM (6% at 300km/s; 2% at 500km/s). 



Appendix C: The derivation of rotation curves 

To fix notations and assumptions we briefly describe how 
rotational velocities arc derived for inclined plane disks. 
The disk may recede with a heliocentric radial velocity 
v£ el . Firstly, coordinates x and y arc defined in the sky 
plane, with their axes along the apparent major and mi- 
nor axis, respectively. The corresponding coordinates in 
the galactic plane are £ = x and i] — y/ cost, where i 
ist the galactic disk's angle of inclination as measured 
from face-on. Galactocentric distances in the sky plane 
are p = \J x 2 + y 2 and correspond to galactocentric dis- 
tances R — + rf — ^/x 2 + (y/ cost) 2 in the plane of 
the galaxy. A polar angle in the plane of the galaxy 
(cos0 = £/R) corresponds to 9 in the sky plane via 
cos0 = cosi/-\/(cosi) 2 + (tan6>) 2 . 

To obtain the observed line-of-sight velocity compo- 
nent v los of a circular (rotational) velocity v rot in the galac- 
tic plane, one first projects v rot with onto a direction 
parallel to the minor axis of the galaxy and then with i 
onto the l.o.s. so that v los = v rot cos sin i or, with the 



above relation between and 9 (and v 



los 



obs 
J hcl 



K 



obs 
1 



V'(cosi) 2 + (tan 0) 2 



sin i cos % 



(C.l) 
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